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ABSTRACT 


The  heat  transfer  along  the  axis  of  a  pulse  detonation  combustor  has  been  characterized 
for  various  frequencies  and  fill  fractions  at  2.5  atmospheres  of  pressure  for  chamber 
refresh  conditions.  In  a  pulse  detonation  combustor,  a  supersonic  detonation  wave  is  the 
method  for  transforming  chemical  energy  into  mechanical  energy  and  the  wave 
propagates  much  faster  than  the  subsonic  flames  in  devices  such  as  rockets  and  ramjets. 
The  flow  field  inside  a  pulse  detonation  combustor  is  highly  turbulent,  unsteady,  and 
varies  largely  during  each  combustion  cycle.  By  determining  the  heat  transfer  properties 
at  multiple  axial  locations  and  the  associated  combustor  wall  temperatures,  proper 
combustor  material  selection  can  ensure  the  material  properties  will  not  deteriorate  and 
therefore  allow  for  practical  operational  lifetimes.  Experimental  testing  measured  the 
axial  heat  transfer  characteristics  in  a  pulse  detonation  combustor  at  various  operating 
conditions  and  multiple  cooling  jacket  locations.  Computer  simulations  were  used  to 
model  the  heat  transfer  inside  the  pulse  detonation  combustor  and  correlate  those 
predications  with  empirical  data.  The  acquired  data  from  the  comparison  of  the  computer 
simulations  and  the  experimental  results  was  correlated  and  demonstrated  good 
agreement.  The  determined  values  should  allow  designers  the  ability  to  consider 
regenerative  fueling  strategies  for  future  systems. 
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I.  INTRODUCTION 


Pulse  detonation  combustion  (PDC)  is  a  technology  in  which  there  has  been  much 
research  and  growing  interest  in  the  last  several  years.  This  technology  is  still  in  the 
developmental  stages,  and  is  demonstrating  the  ability  to  produce  improved  performance 
over  many  proven  sources  of  propulsion  and  power  generation.  Technical  challenges 
associated  with  PDC  systems  are  being  overcome  and  are  revealing  the  potential 
advantages  over  current  technologies.  If  applied  to  flight,  pulse  detonation  combustion 
can  be  used  to  power  a  tactical  missile  with  a  Pulse  Detonation  Engine  (PDE)  and 
provide  greater  range  for  the  same  amount  of  fuel  as  compared  to  other  current 
technologies. 

As  advances  continue  in  material  sciences,  the  PDE  could  see  a  wider  operational 
range  at  higher  flight  Mach  numbers  by  using  lighter  and  stronger  materials  that  could 
withstand  higher  operating  temperatures  and  conditions  within  the  combustor.  With 
significant  improvements  in  thennal  efficiency  and  simplicity  of  design  as  compared  to 
other  types  of  engines,  PDEs  could  make  an  excellent  propulsion  source  for  supersonic 
tactical  missile  designs. 

As  seen  in  Figure  1,  the  specific  impulse,  I  ,  or  fuel  efficiency,  of  a  PDE  is 

greater  than  ramjets,  scramjets,  ducted  rockets  and  solid  rockets.  The  upper  band  of  a 
turbojet  represents  fuel-lean  operation  and  exceeds  a  PDE  in  specific  impulse.  The  PDE 
performance  could,  therefore,  be  improved  substantially  if  operated  fuel-lean  or  under 
partial  fill  conditions.  A  PDE,  however,  is  constructed  with  very  few  moving  parts  as 
compared  to  a  turbojet,  therefore  greatly  reducing  costs  associated  with  production  and 
construction. 
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Figure  1 .  Comparison  of  High-Speed  Propulsion  Technologies  (From  [1]) 


Further  research  is  also  ongoing  in  using  hybrid  PDC  gas  turbine  systems  for 
power  generation  applications.  Gas  turbines  are  known  for  great  overall  efficiencies  but 
often  only  while  operating  at  design  conditions,  which  are  where  most  U.S.  naval  ships 
do  not  spend  extensive  time  operating.  Gas  turbines  lose  their  great  efficiency  when 
operating  at  lower  speeds,  while  ships  are  patrolling  or  while  in  port  at  idle  settings.  By 
using  a  hybrid  PDE  gas  turbine  system,  an  efficient  PDC  could  improve  performance  at 
lower  rotational  speeds,  therefore  increasing  the  overall  system  efficiency  during  low  and 
high  power  demand  operations. 

In  order  to  achieve  improved  performance  in  detonation-based  systems,  a  full 
understanding  of  the  thermal  loads  experienced  by  a  PDC  is  essential  for  developing  their 
practical  use.  In  a  PDC,  the  combustor  is  operated  in  a  cyclical  manner  by  producing 
multiple  detonation  waves  per  second.  These  detonation  waves  create  substantial 
transient  chamber  pressures  at  high  operating  temperatures.  The  extremely  transient  and 
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brutal  environment  makes  it  especially  difficult  to  experimentally  obtain  time  resolved 
heat  transfer  measurements  and  few  studies  have  been  published  in  this  area. 

This  becomes  a  challenging  problem  to  solve,  which  requires  comprehension  of 
fluid  dynamics  and  heat  transfer  of  turbulent  flow.  The  turbulence  created  in  the 
combustor  creates  random  fluctuations  in  the  fluid  and  on  both  a  global  and  a  local  scale 
the  flow  becomes  highly  unsteady.  Hence,  a  time  averaged  value  may  be  used  in  order  to 
achieve  a  parameter  that  is  independent  of  time  and  assumed  to  be  quasi  steady.  Naples, 
Hoke,  and  Schauer  [2]  at  the  Air  Force  Research  Laboratory  performed  an  experiment 
that  studied  the  heat  loads  from  steady  deflagration  and  pulsed  detonation  combustion  by 
placing  a  cross  flow  tube  within  the  flow  field  of  the  combustor  in  order  to  measure  the 
time  averaged  heat  transfer.  The  unsteady  flow  of  a  PDE  cycle  develops  an  environment 
that  is  thermally  different  from  that  of  a  continuous  flow  cycle.  The  results  displayed  a 
20-30%  lower  heat  load  in  pulse  detonation  as  compared  to  steady  deflagration.  The 
experiment  also  revealed  that  the  un-insulated  combustor  wall  temperature  for  steady 
state  deflagration  was  in  excess  of  1,400°F  while  the  pulse  detonation  cycle  resulted  in  an 
approximate  temperature  of  1,050°F.  Concluding,  the  violent  combustion  of  pulse 
detonation  transferred  less  heat  to  the  wetted  surface  area  than  a  steady  deflagration 
combustor  setup. 

At  the  NASA  Glenn  Research  Center,  Paxson  [3]  conducted  research  by 
comparing  experimentally  measured  and  numerically  simulated,  time-averaged,  point 
heat  transfer  rates  in  a  PDE  using  a  cylinder  in  cross-flow  and  a  spool  design.  The  cross- 
flow  cylinder  was  placed  directly  in  the  combustion  gas  flow  field,  which  produces 
different  results  as  compared  to  the  heat  transfer  occurring  through  the  combustor  wall. 
The  computational  fluid  dynamics  (CFD)  model  was  validated  by  agreement  with 
measured  and  computed  results  of  several  locations  under  different  operating  conditions. 
These  results  were  completed  on  an  open-ended  combustion  chamber  that  allowed  flow 
reversals  from  the  atmosphere  to  enter  the  combustion  chamber,  which  provided  cooling 
that  could  not  be  accurately  accounted  for. 

Further  research  is  needed  to  characterize  the  changes  in  heat  transfer  rates  along 

the  length  of  the  combustor  without  introducing  any  obstacles  into  the  already 
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complicated  flow  field  and  with  a  nozzle  ended  combustion  chamber  to  prevent  flow 
reversals.  These  results  would  lead  to  time-averaged  values  for  an  effective  heat  transfer 
coefficient  of  the  combustor.  Ultimately,  it  is  desirable  to  obtain  local  parameters  at  the 
combustor  wall  vs.  time-averaged  values. 

The  objective  of  this  research  is  to  characterize  the  heat  transfer  properties  along 
the  axis  of  a  PDC  for  various  frequencies  and  fill  fractions.  Experimental  testing  was 
performed  to  determine  the  heat  transfer  characteristics  at  multiple  axial  combustor 
locations  for  operations  at  20,  30  and  40  Hz  with  varying  fill  fractions  of  14,  16,  and  18 
milliseconds.  Experimental  testing  contributed  in  detennining  the  regions  of  highest  heat 
transfer  and  hence  where  the  greatest  amount  of  cooling  is  required.  Computational  fluid 
dynamics  modeling  was  used  to  model  the  heat  transfer  inside  a  combustor  during  a  pulse 
detonation  event,  and  correlated  with  empirical  data. 
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II.  BACKGROUND 


A.  DESCRIPTION  OF  BASIC  PULSE  DETONATION  ENGINE  OPERATION 

A  prototypical,  straight  -tube,  valve-less  PDC  is  closed  at  one  end  and  opens  at 
the  other  end.  In  some  cases,  a  nozzle  may  or  may  not  be  included  as  part  of  the  PDC. 
The  various  stages  of  the  pulse  detonation  combustion  cycle  consist  of  four  basic  steps: 
detonation,  blow  down,  purge  and  fill.  Figure  2  displays  the  basic  combustion  cycle  of  a 
PDE  and  give  a  visual  understanding  of  the  events  that  are  occurring  during  a  cycle. 


Simple  basic  cycle,  implementation  technically  very  challenging 


Figure  2.  Simplified  Ideal  PDE  Operation  Cycle  (After  [4]) 


In  the  first  step  during  the  fill  stage,  the  cycle  is  prepared  by  injecting  a  mixture  of 
fuel  and  oxidizer  into  the  head  end  of  the  combustor.  In  step  two,  the  mixture  is  ignited, 
which  creates  a  deflagration  event  in  the  combustion  chamber.  Steps  three  and  four 
displays  the  resulting  flame  that  travels  some  distance  and  then  experiences  deflagration- 
to-detonation  transition  (DDT).  The  detonation  wave  travels  down  the  combustor  and 
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exits  to  the  atmosphere.  Blow  down  starts  to  occur  in  step  five  after  the  detonation  wave 
exits  the  tube  and  expansion  waves  move  upstream  towards  the  end-wall.  Once  these 
expansion  waves  reach  the  end-wall,  the  pressure  in  the  chamber  begins  to  drop  while 
still  producing  a  decreasing  thrust.  In  the  final  step,  the  end-wall  pressure  has  dropped 
low  enough  to  allow  purge  gases  to  flow  into  the  combustion  chamber.  During  the  purge 
stage,  a  slug  of  purge  gas  is  typically  injected  to  isolate  the  combustion  products  from  the 
next  fill  and  then  a  new  cycle  begins. 

B.  COMBUSTION  PROCESSES 

1.  Deflagration 

The  most  common  type  of  combustion  process  is  the  deflagration  combustion 
process  that  exists  in  majority  of  engines  in  use  today.  The  majority  of  all  flight  engines 
and  ground-based  power  generation  systems  bum  fuel  as  a  deflagration  to  release  the 
energy  contained  within  the  fuel.  During  the  deflagration  process,  combustion  is 
occurring  at  subsonic  velocities.  This  process  occurs  at  nearly  constant-pressure 
conditions.  The  combustion  wave  propagates  down  the  combustor  and  energy  is 
transferred  to  the  working  fluid  via  thermal  diffusion.  Deflagration  combustion  can  be  a 
continuous  process  such  a  burning  open  flame. 

2.  Detonation 

The  combustion  process  of  a  PDE  transitions  from  a  deflagration  combustion 
process  to  detonation  combustion.  Understanding  the  detonation  process  is  essential  in 
evaluating  the  heat  transfer  process  of  the  combustion  chamber.  Detonation  occurs  at 
supersonic  combustion  wave  speeds  during  a  constant  volume  process.  A  fuel/air 
mixture  traveling  at  a  subsonic  velocity  is  compressed  by  a  normal  shock  wave 
propagating  at  a  high  Mach  number,  which  is  then  followed  by  a  rapid  release  of  heat  and 
a  sudden  rise  in  pressure.  The  coupling  of  a  strong  shock  to  a  region  of  heat  release  is 
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called  a  detonation.  This  is  a  very  violent,  rapid,  exothermic  reaction,  which  produces  a 
harsh  environment  in  which  it  is  difficult  to  experimentally  document  the  instantaneous 
changes  in  local  pressures  and  temperatures. 

The  detonation  process  results  in  pressure  gains  across  the  combustion  wave 
where  a  deflagration  combustion  process  results  in  a  minimal  or  no  pressure  rise.  Also 
through  detonation,  combustion  products  exhibit  higher  temperatures  due  to  the  inherent 
compression  and  less  dissociation.  A  deflagration  will  often  occur  as  a  lower  flame 
temperature  and  the  process  is  typically  continuous.  During  detonation,  the  higher  flame 
temperatures  are  only  periodic  from  the  nature  of  the  cycle.  The  fill  and  purge  steps 
result  in  a  cooling  portion  of  the  cycle. 

Detonation  of  a  fuel/air  mixture  is  often  difficult  to  initiate  within  a  short  length 
combustor.  One  method  to  achieve  detonation  is  to  start  with  deflagration  combustion 
and  then  transition  the  flame  front  to  a  detonation  to  by  inserting  obstacles  in  the  flow 
field.  The  obstacle  causes  the  flow  to  speed  up  as  well  as  creating  turbulent  mixing.  This 
results  in  a  process  known  as  Deflagration  to  Detonation  Transition  (DDT).  A  common 
obstacle  used  in  the  flow  field  is  a  Shchelkin  spiral,  but  this  causes  a  noticeable  pressure 
loss  during  the  filling  and  detonation  portion  of  the  cycle  [5].  A  more  recent  study  by 
Dvorak  [6]  demonstrates  that  using  obstacles  such  as  swept  ramps  can  provide  at  least  a 
27%  improvement  over  the  total  pressure  loss  of  a  wall  spiral  with  the  same  DDT 
performance. 

C.  HEAT  TRANSFER  OF  THE  COMBUSTION  CHAMBER 

Extended  operation  of  a  PDC  at  practical  frequencies  and  pressures  requires 
cooling  of  the  combustor  wall.  A  water  cooling  system  is  currently  in  use  at  the  Naval 
Postgraduate  School  Rocket  Propulsion  Laboratory  to  ensure  overheating  and  failure  of 
the  combustor  wall  does  not  occur  and  allows  overall  heat  transfer  values  to  be 
detennined  for  future  design  purposes.  Heat  is  transferred  from  the  bulk  of  the 
combustion  gases  in  the  combustion  chamber  to  the  combustor  wall  surface  via  forced 
convection.  From  there,  the  heat  flows  through  the  combustor  wall  to  inner  wall  of  the 
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cooling  jacket  by  means  of  conduction.  Finally,  forced  convection  occurs  as  the  cooling 
water  bulk  fluid  transports  heat  away  from  the  surface  of  the  cooling  jacket  inner  wall. 

1.  Conduction  Through  the  Combustor/Cooling  Jacket  Wall 

Heat  is  transferred  by  conduction  through  the  combustor  wall  under  the  assumed 
one-dimensional,  steady  state  conditions.  Heat  transfer  is  occurring  all  the  way  through 
the  wall  without  any  internal  generation  of  thermal  energy.  As  shown  in  Figure  3,  heat  is 
transferred  in  the  x  direction  and  temperature  is  a  function  of  x.  The  combustor  wall  is 
physically  separating  the  hot  combustor  gases  from  the  cold  cooling  water.  Heat  transfer 
occurs  from  the  hot  gas  through  the  wall  to  the  cold  liquid.  If  both  the  gas  wall  and  liquid 
wall  temperatures  are  known,  the  conduction  heat  transfer  rate  can  be  detennined,  as  well 
as  a  temperature  distribution.  The  inherently  unsteady  conditions  of  the  detonation 
process  is  the  driving  force  for  heat  transfer  to  occur  across  the  combustor  wall  to  a  quasi 
steady  state  process  of  convection  provided  by  the  cooling  water  side  of  the  combustor 
wall. 
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(1) 


—(A'—)  =  0 
dx  dx 


For  a  one-dimensional,  steady  state  conduction  in  a  plane  wall  with  no  heat 
generation  and  a  constant  thermal  conductivity,  the  temperature  through  the  wall  will 
vary  linearly  as  x  increases  [5].  The  conduction  heat  transfer  rate,c/v,  and  heat  flux  ,q"x , 
through  a  plane  wall  is  given  by  Equation  (2)  and  (3). 
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The  wall  area  nonnal  to  the  direction  of  heat  transfer  is  noted  as  the  cross- 
sectional  area,  A.  The  combustor  wall  is  made  for  4340  steel  that  has  a  thermal 
conductivity,  k,  of  44.5  W/m-K  at  300  K.  The  temperature  of  the  combustor  inner  wall 
T  waii ,  is  extremely  difficult  to  experimentally  measure  due  to  the  high  pressure, 

temperatures  and  transient  environment.  The  small  clearances  in  the  cooling  jacket  make 
it  difficult  to  measure  the  liquid  side  of  the  inner  jacket  wall  temperature  Tlwall,  so  an 

energy  balance  with  the  cooling  water  is  required  in  order  to  determine  the  inner  wall 
cooling  jacket  temperature. 


2.  Convection  Along  the  Cooling  Jacket  Inner  Wall 

The  primary  heat  transfer  between  a  surface  and  a  fluid  moving  over  the  surface  is 
known  as  convection.  For  a  transfer  of  energy  to  occur,  the  fluid  that  may  exist  as  a 
liquid  or  a  gas  moves  at  a  prescribed  velocity  and  there  exists  a  temperature  difference 
between  the  temperature  of  the  surface  and  the  temperature  of  the  moving  fluid.  The 
local  heat  flux,  q” ,  is  expressed  in  Equation  (4). 

<?'  =  -A  =KT,.„  -TJ  (4) 
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The  local  convection  coefficient,  h,  will  vary  on  a  surface  as  the  flow  conditions 
vary  from  point  to  point.  This  variation  in  the  local  convection  coefficient  will  cause  the 
local  heat  to  fluctuate  across  the  surface  as  well.  By  using  an  average  convection 
coefficient,  h  ,  for  the  whole  surface,  then  the  total  heat  transfer  rate,  q,  can  be  expressed 
by  Equation  (5). 

q  =  hAs(T-TJ  (5) 

Determination  of  the  local  and  average  convection  coefficients  is  not  an  easy 
problem.  The  coefficients  depend  on  several  fluid  properties  such  as  density,  viscosity, 
thennal  conductivity  and  specific  heat,  as  well  as  the  flow  conditions  and  surface 
geometries  [7]. 


a.  Velocity  Boundary  Layer 

Figure  4  illustrates  the  flow  over  a  flat  plate.  The  wall  surface  where  the 
fluid  particles  make  contact,  the  fluid  velocity  is  assumed  to  be  zero  and  is  often  referred 
to  as  a  no-slip  condition.  These  near-zero  velocity  particles  hinder  the  fluid  layer  that  is 
adjacent  to  the  wall.  This  hindrance  between  layers  continues  some  distance  away  from 
the  wall  resulting  in  layers  of  shear  stresses  until  the  distance  from  the  wall  increases  to  a 
value  where  the  velocity  of  the  fluid  changes  to  that  of  the  free  stream.  To  properly  solve 
a  convection  problem,  it  must  be  determined  whether  the  boundary  layer  is  laminar  or 
turbulent.  The  laminar  boundary  layer  fluid  motion  is  ordered  and  in  the  same  direction. 
The  turbulent  boundary  layer  grows  larger  than  the  laminar  boundary  layer  and  is  erratic 
with  fluctuations  in  all  directions. 
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Figure  4.  Velocity  Boundary  Layer  Development  on  a  Flat  Plate  (From  [7]) 


As  the  boundary  transitions  from  laminar  to  turbulent,  the  boundary  layer 
grows  larger,  while  the  shear  stress  and  local  heat  transfer  coefficient  significantly 
increase.  Notice  in  Figure  5  the  sudden  increase  in  the  local  heat  transfer  coefficient  as 
flow  transitions  from  laminar  to  turbulent. 


Transition 

Figure  5.  Variation  of  Velocity  Boundary  Layer  Thickness  5  and  the  Local  Heat 
Transfer  Coefficient  h  for  Flow  Over  an  Isothennal  Flat  Plate  (From  [7]) 
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b.  Thermal  Boundary  Layer 


Similar  to  the  growth  of  a  velocity  boundary  layer,  a  thermal  boundary 
develops  as  long  as  there  is  a  difference  in  the  plate  surface  temperature  and  the  free 
stream  fluid  temperature.  Thennal  equilibrium  is  achieved  at  the  wall  surface  between 
the  wall  and  the  zero  velocity  particles.  At  this  location,  energy  transfer  only  occurs  by 
conduction  since  there  is  zero  fluid  motion.  Energy  is  then  transferred  by  convection 
from  the  stationary  particles  to  the  adjacent  particle  in  the  moving  fluid  layer  and 
continues  to  transfer  through  layers.  The  growth  of  the  thermal  boundary  layer  is 
illustrated  in  Figure  6.  The  rate  of  heat  transfer  is  determined  by  the  conditions  in  the 
thermal  boundary  layer  and  influences  the  heat  flux  and  the  local  convective  heat  transfer 
coefficient. 


Figure  6.  Thermal  Boundary  Fayer  Development  on  an  Isothermal  Flat  Plate 

(From  [7]) 


c.  Fully  Developed  Internal  Flow 

In  internal  flow,  the  fluid  is  confined  by  a  surface;  therefore,  boundary 
layer  development  is  limited  due  to  this  constraint.  The  point  at  which  the  velocity 
profde  no  longer  changes  is  then  considered  fully  developed.  For  flow  through  the 
cooling  jacket  channels,  it  must  be  determined  whether  the  flow  is  laminar  or  turbulent  in 
order  to  establish  at  which  point  flow  is  fully  developed.  For  internal  flow,  the 
commencement  of  turbulence  begins  at  a  Reynolds  number  of  approximately  2300,  while 
numbers  as  large  as  10,000  are  needed  in  order  to  reach  fully  turbulent  conditions.  For  a 

circular  tube,  the  Reynolds  number  can  be  determined  from  Equation  (6). 

12 


(6) 


ReD 


M 


Inviscid  flow  region  „ 

6  Boundary  layer  region 

«(r.  x) 


Figure  7.  Laminar,  Hydrodynamic  Boundary  Layer  Development  in  a  Circular  Tube 

(From  [7]) 


Since  there  is  a  variation  in  the  velocity  through  the  cross  section  of  the 
pipe,  the  mean  velocity,  um,  is  used  to  solve  problem  with  internal  flows.  Mass  flow  rate 
can  then  be  determined  from  the  mean  velocity,  density  and  cross-sectional  area  as  in 
Equation  (7). 

m  =  pumAc  (7) 


For  a  circular  tube  with  incompressible  flow  with  a  constant  mass  flow 
rate,  the  Reynolds  number  is  equivalent  to  Equation  (8).  The  Reynolds  number  is 
interpreted  as  the  ratio  of  the  inertia  and  viscous  forces. 


Re0 


4m 

nD/u 


(8) 


Since  the  cooling  jacket  water  flow  channels  are  square,  noncircular  tubes 
use  an  effective  diameter  known  as  the  hydraulic  diameter,  shown  in  Equation  (9).  For 
the  cooling  jackets,  the  hydraulic  diameter  is  to  be  used  in  calculations  for  Reynolds  and 
Nusselt  numbers.  In  Equation  (9),  P  is  the  wetted  perimeter. 


Du 


(9) 
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The  Prandtl  number  is  required  to  calculate  the  Nusselt  number.  The 
Prandtl  number  is  interpreted  as  the  ratio  of  momentum  and  thermal  diffusivities. 


Pr  = 


(10) 


The  Nusselt  number  can  be  determined  from  the  Dittus-Boelter  equation 
with  the  limitations  listed  below.  The  Nusselt  number  is  a  dimensionless  temperature 
gradient  of  the  surface. 


Nud  =  0.023  Re^  Pr" 

0.7  <  Pr  <  160 

Re^  >  10,000 

—  >10 

ID  J 

(Ts  >  Tm)  =>  n  =  0.4 


(11)  (From  [7]) 


The  average  heat  transfer  coefficient,  h  ,  can  be  detennined  from  the 
average  Nusselt  number. 

Nu  =  —  (12) 

kf 

Similar  to  Equation  (4),  the  mean  temperature,  Tm,  is  used  as  a  reference 
temperature  for  internal  comparable  to  the  free  stream  temperature  for  external  flows. 
Newton’s  law  of  cooling  can  be  expressed  as  in  Equation  (13). 

q"  =  h(Ts  —Tm)  (13) 

The  total  cooling  jacket  heat  transfer  rate  can  be  determined  by  using  the 
specific  heat  of  the  cooling  water  and  the  change  in  temperature  from  the  inlet  to  outlet  of 
the  cooling  jacket.  Equation  (14)  is  an  appropriate  general  equation  that  applies 
regardless  of  the  nature  of  the  surface  thermal  or  tube  flow  conditions  [7]. 

Vconv  =  ™Cp(Tm,ou<  -  T,n,J  (14) 
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The  mean  temperature  in  Tm  in ,  and  the  mean  temperature  out  Tm  out ,  are 

measured  by  a  thermocouple.  The  mass  flow  is  physically  measured;  therefore,  the  heat 
flux  is  determined  by  dividing  the  heat  transfer  rate  by  the  inner  surface  area  of  the 
cooling  jacket.  After  the  Nusselt  number  is  determined,  the  average  heat  transfer 
coefficient  is  found  by  using  Equation  (12).  The  surface  temperature  of  the  cooling 
jacket  wall  is  then  evaluated  using  the  heat  flux  from  Equation  (14).  By  knowing  the 
surface  temperature  of  the  cooling  jacket  wall,  the  surface  temperature  on  the  inside  wall 
of  the  combustor  is  resolved  by  conduction  in  Equation  (3). 

3.  Convection  Along  the  Combustor  Inner  Wall 

Convection  along  the  inner  wall  of  the  combustion  chamber  is  very  difficult  to 
evaluate  based  on  the  continuously  changing  flow  conditions  due  to  detonation  and 
refresh  periods.  Even  if  the  combustor  always  had  a  steady  flow,  in  order  to  consider 
fully  developed  turbulent  flow,  the  length  of  the  combustor  would  have  to  be  ten  times 
the  hydraulic  diameter.  Since  the  diameter  of  the  combustion  chamber  is  7.62  cm,  a 
L/D  ratio  greater  than  ten  would  require  the  combustor  length  to  be  greater  than  30”  in 
length.  This  condition  is  not  met,  therefore,  the  Nusselt  number  calculation  in  Equation 
1 1  cannot  be  used  to  aid  in  calculating  the  heat  transfer  coefficient.  The  thermal 
boundary  layer  thickness  is  zero  at  the  tube  entrance  and  the  convection  coefficient 
would  be  extremely  large  at  this  location.  As  the  thermal  boundary  develops,  the 
convection  coefficient  decays  rapidly. 
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Figure  8.  Thermal  Boundary  Layer  Development  in  a  Heated  Circular  Tube 

(From  [7]) 


Since  the  heat  flux  across  the  tube  wall  is  known,  an  alternate  approach  using 
Equation  13  can  be  used  to  determine  the  effective  heat  transfer  coefficient  of  the  inside 
combustor  wall.  One  expectation  is  that  since  detonation  does  not  occur  until  after  the 
installed  ramps,  the  observed  heat  transfer  rate  will  vary  axially  depending  upon  the 
location  inside  the  combustion  chamber. 
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III.  DESIGN  AND  EXPERIMENTAL  SETUP 


A.  PULSE  DETONATION  ENGINE 


The  experimental  testing  occurred  at  the  Rocket  Propulsion  Laboratory  (RPL),  an 
off-campus  testing  facility  that  is  part  of  the  Naval  Postgraduate  School  in  Monterey, 
California.  The  PDC  and  a  majority  of  the  equipment  used  for  testing  were  located  in 
Test  Cell  #2.  Additional  equipment  was  located  in  Test  Cell  #3  and  in  the  Control  Room. 
The  test  equipment  included  the  combustor  section,  cooling  jackets,  ethylene/air  delivery 
systems,  cooling  water  system,  data  acquisition  system,  and  PDC  controller.  Most  of  the 
equipment  was  used  in  previous  experimentations;  however  additional  thermocouples 
and  pressure  transducers  were  added  to  accurately  monitor  various  water  jacket 
conditions. 


Figure  9.  Naval  Postgraduate  School  Rocket  Propulsion  Laboratory  Test  Cell  #2 


The  valve-less  single  tube  PDC  design  at  NPS  has  continuously  evolved  through 
research  by  graduate  students  for  the  past  six  years.  The  PDC  consists  of  a  fuel  and  air 
injection  system,  an  ignition  system,  and  a  combustion  chamber  that  is  coupled  with  a 
cooling  water  system.  The  fuel  injector  and  ignition  system  is  connected  to  the  PDC 
controller  and  all  of  the  diagnostics  are  recorded  by  a  data  acquisition  system  to 
determine  performance. 
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B.  FUEL  AND  AIR  DELIVERY 

1.  Air 

The  supply  air  was  delivered  to  the  PDC  by  a  5.08  cm  diameter  supply  line  after 
being  heated  by  a  vitiator.  This  allowed  for  simulation  of  various  combustor  inlet 
conditions.  The  vitiator  operated  by  injecting  hydrogen  into  the  main  air  supply  flow, 
which  was  then  ignited  with  a  hydrogen/air  torch.  Hot  air  from  the  vitiator  heated  the 
PDE  downstream  piping  and  combustor  to  the  desired  conditions.  This  was 
accomplished  by  running  the  vitiator  for  approximately  15  seconds.  Air  temperature 
remained  steady  for  several  seconds  after  the  vitiator  was  turned  off,  but  cooler  incoming 
air  would  be  experienced  by  the  PDE  for  longer  duration  runs. 

The  heated  air  was  split  and  delivered  to  the  combustor  through  four  3.81  cm 
fueling  arms  where  fuel  was  added  to  the  air.  In  each  of  the  inlet  arms,  chokes  were 
installed  in  order  to  separate  the  vitiator  from  downstream  pressure  oscillations.  By 
splitting  the  air  and  fuel  into  four  separate  arms,  proper  fuel/air  mixing  occurred  while 
providing  a  more  even  injection  of  the  fuel/air  mixture  into  the  combustor. 

2.  Ethylene 

Controlled  delivery  of  the  ethylene  fuel  in  a  time -varying  manner  is  necessary  in 
order  to  supply  proper  stoichiometry  of  the  fuel/air  mixture  provided  to  the  combustor. 
Varying  the  mass  flow  rate  of  the  injected  fuel  was  accomplished  by  varying  the  supply 
pressure  with  Tescom  regulators.  The  varied  fuel  injection  pressure  adjusts  the  mass 
flow  rate  of  the  fuel  thereby  dictating  the  equivalence  ratio  (cp)  of  the  mixture  for  each 
engine  cycle.  In  Equation  15,  the  equivalence  ratio  is  determined  by  dividing  the  mass 
ratio  fuel  to  air  used  in  testing  by  the  mass  ratio  of  fuel  to  air  for  the  stoichiometric 
condition. 
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The  highest  thrust  values  occur  when  the  equivalence  ratio  of  the  fuel/air  mixture 
is  near  unity.  If  the  mixture  is  rich  or  an  equivalence  ratio  greater  than  unity,  fuel  will 
left  unburned  since  there  is  more  fuel  that  can  be  combusted  for  the  given  amount  of 
oxidizer.  Lower  specific  impulse  is  the  result  of  the  wasted  fuel  while  still  producing 
similar  thrust  levels.  In  opposition,  a  lean  mixture  that  yields  an  equivalence  ratio  less 
than  unity  result  in  lower  thrust  values,  but  can  generate  higher  fuel-based  specific 
impulses. 

(16) 

(17) 


Ethylene  is  delivered  to  the  PDC  from  an  accumulator  installed  prior  to  the 
research  performed  by  Nichols  [8].  The  installed  accumulator  allows  for  constant 
pressure  delivery  of  ethylene  and  longer  duration  runs  of  the  PDC  in  which  both  are 
required  to  generate  steady-state  heat  transfer  condition.  The  blue  accumulator  in  Figure 
10  has  a  piston  that  is  pressured  inside  a  cylindrical  pressure  vessel.  Ethylene  fills  the 
accumulator  on  one  side  of  the  piston,  and  nitrogen  gas  pressurizes  the  ethylene  by 


SP  (™Pg0 )  (mg0) 


Figure  10.  Test  Cell  2  Fuel  and  Gas  Tanks 
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supplying  the  required  pressure  on  the  other  side.  Pressure  of  the  nitrogen  gas  was 
controlled  by  Tescom  regulators.  The  pressurized  ethylene  is  supplied  to  the  PDC  at  the 
required  pressure. 


Four  independent  injectors  located  in  four  fuel  arms  were  used  to  supply  ethylene 
to  the  PDE.  A  common  feed  manifold  with  four  electronically-controlled  high  frequency 
Valvetech  (PN#  15060-2)  solenoid  valve  injectors  were  mounted  to  the  fuel  arms 
downstream  of  the  flow  chokes.  Supplied  air  was  mixed  with  the  gaseous  fuel  mixture 
prior  to  entering  the  combustion  chamber.  One  of  the  two  new  fast  response  Vavletech 
(PN#12177-2)  solenoid  valve  injectors  used  by  Nichols  [8]  was  defective  and  returned  to 
the  factory  for  service. 


Figure  1 1 .  PDE  Fuel  Arms  and  Injectors 


3.  Ignition  System 

A  small-scale  Transient  Plasma  Ignition  (TPI)  system  was  used  to  control  the 
ignition  of  the  PDC.  A  TPI  box  generates  a  high  voltage  signal  to  an  electrode  assembly. 
The  electrode  assembly  is  inserted  into  the  combustion  chamber  through  a  machined 
orifice.  The  BNC  500  Pulse  generator  pulsates  at  a  preset  frequency  to  the  BNC  575  in 
order  to  produce  two  output  waveforms  for  each  pulse.  These  waveforms  provide  a 
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signal  to  produce  a  rapid  charge  input  to  the  High  Volt  Pulse  Generator.  The  electrode 
then  develops  a  corona  discharge  to  the  electrode  shroud.  The  TPI  cycle  is  displayed  in 
Figure  12. 


Figure  12.  Flow  Path  for  Transient  Plasma  Ignition  (From  [6]) 


The  electrode  shroud  pictured  on  the  left  in  Figure  13  was  used  on  multiple  runs 
of  the  PDC.  During  runs  of  five  seconds  or  longer,  flame  holding  was  obviously  present 
in  the  PDC  exhaust.  The  long  shroud  was  suspect  of  causing  auto  ignition  of  the  fuel/air 
mixture  due  to  the  high  shroud  temperatures  during  extended  runs.  The  shroud  was  cut  to 
various  shorter  lengths  to  minimize  the  potential  for  flame  holding  to  occur  and  did 
improve  the  engine  operation. 


Figure  13.  TPI  Electrode  Shroud 
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c. 


COMBUSTOR  ASSEMBLY 


1.  Combustor 

The  combustor  assembly  consisted  of  a  7.62  cm  diameter  combustion  tube  with  a 
convergent  divergent  nozzle  on  the  end.  The  eight  combustion  tube  sections  and  the 
nozzle  were  assembled  together  to  form  the  complete  combustor  assembly.  Each  section 
was  made  of  4340  steel  and  has  a  cooling  jacket  that  is  supplied  by  the  cooling  water 
system. 

The  primary  combustion  tube  consisted  of  eight  7.62  cm  sections  for  an  overall 
length  of  0.6096  m.  The  combustor  section  had  a  7.62  cm  inside  diameter  and  it  was 
constructed  from  4340  annealed  steel.  Each  section  was  a  two  piece  assembly  designed 
by  Thomas  Lipoma,  a  NPS  RPL  summer  intern  student.  The  inner  tube  was  pressed  into 
an  outer  casing  with  flanges  on  both  ends  to  allow  each  section  to  be  bolted  together.  An 
O-ring  was  inserted  into  the  groove  on  the  aft  flange  face  to  insure  proper  sealing 
between  each  section.  The  combustor  had  an  inner  wall  thickness  of  0.3175  cm  to 
separate  combustion  gases  from  the  cooling  water. 


Figure  14.  PDE  Combustor 
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Each  section  had  a  hole  drilled  in  the  aft  flange  that  penetrated  through  the 
combustor  wall.  The  holes  were  180  degrees  apart  at  the  top  and  bottom  of  the 
combustor  tube.  Each  hole  allowed  for  the  installation  of  an  obstacle  onto  the  combustor 
inner  wall.  The  obstacles  are  discretely  placed  along  the  combustor  wall  to  achieve  DDT. 
The  preferred  obstacle  shape  to  achieve  a  balance  between  detonation  initiation 
performance  and  pressure  loss  was  determined  by  research  conducted  by  Dvorak  [6]. 
The  obstacle  designed  was  a  swept  ramp  in  a  configuration  known  as  2R.  180.4S  [6]  as 
seen  in  Figure  15. 


SWEPT  TALL 

Figure  15.  Dvorak  Swept-Tall  Ramp  (From  [6]) 

This  research  used  four  sets  of  ramps  place  in  combustor  sections  four  through 
seven.  Each  ramp  was  bolted  in  place  with  a  thin  layer  of  Room  Temperature 
Vulcanizing  (RTV)  silicone  rubber  between  the  ramp  and  the  combustor  wall.  In  the 
sections  that  did  not  utilize  the  use  of  an  obstacle,  the  drilled  holes  were  plugged  with  a 
cap  and  lock  wired  into  place. 


Figure  16.  Swept  Ramp  Configuration  Inside  Combustor  [From  (8)] 
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2.  Cooling  Jackets 


The  channels  in  between  the  inner  tube  and  the  outer  casing  allowed  for  the 
passage  of  cooling  water  and  can  be  seen  in  Figure  17.  A  cooled  combustion  chamber 
was  required  due  to  the  extreme  temperatures  and  pressures  experienced  during  the 
operation  of  the  PDE.  By  cooling  the  combustor,  the  wall  could  maintain  its  integrity 
during  extended  durations.  Part  of  this  research  was  to  determine  the  amount  of  cooling 
required  of  the  combustor  wall  to  achieve  sustained  operation  of  the  PDE. 


Figure  17.  PDE  Cooling  Jacket 


Each  water  cooling  channel  was  0.635  cm  wide  by  0.635  cm  high.  In  each  7.62 
cm  combustor  section,  the  channel  wrapped  around  the  combustor  eight  times  in  order  to 
achieve  a  channel  length  slightly  greater  than  0.762  m.  Each  of  the  eight  sections  was 
attached  to  a  common  supply  manifold  and  a  common  discharge  manifold.  Mass  flow 
rates  were  measured  at  various  differential  pressures  through  each  combustor  section  to 
accurately  determine  the  amount  of  heat  transfer  occurring  through  the  combustor  wall  to 
the  cooling  water. 
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3. 


Cooled  Nozzle 


A  two-piece  nozzle  assembly  similar  to  the  combustor  section  utilized  inner 
channels  to  allow  a  passageway  for  cooling  water.  The  convergent-divergent  nozzle  was 
design  with  a  15  degree  half  angle  and  was  9.21  cm  in  length.  The  nozzles  was  designed 
by  NPS  RPL  intern  student  and  used  in  the  research  conducted  by  Nichols  [8].  A 
detailed  drawing  can  be  viewed  in  Appendix  B. 

4.  Cooling  Water  System 

The  cooling  water  system  provided  cooling  water  through  each  combustor  and  the 
nozzle  section  to  limit  the  combustor  bulk  wall  temperature.  Cooling  water  was 
delivered  to  the  combustor  from  a  115  gallon  cooling  water  tank  via  a  MTH  brand 
regenerative  turbine  pump,  Model  284K  BF  with  a  capacity  of  60  GPM.  The  supply  line 
pressure  maintained  69000  Pa  with  a  20700  Pa  differential  pressure  across  the  supply  and 
discharge  manifolds.  The  position  of  the  recirculation  line  ball  valve  determined  the 
system  operating  pressure  as  well  as  the  differential  pressure  across  the  manifolds.  A 
cooling  water  system  diagram  is  included  in  Appendix  B. 


Figure  18.  Cooling  Water  Tank  (Left)  and  Regenerative  Turbine  Water  Pump  (Right) 

(From  [8]) 
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The  cooling  water  system  delivered  the  required  mass  flow  rate  in  order  to 
achieve  proper  cooling  of  the  exterior  combustor  wall  by  means  of  forced  convection.  A 
common  thermocouple  and  pressure  transducer  were  installed  on  the  supply  manifold. 
On  the  discharge  manifold,  individual  thermocouples  were  installed  to  monitor  the  outlet 
temperature  of  each  cooled  combustor  section  and  the  cooled  nozzle.  A  differential 
pressure  cell  was  connected  between  both  manifolds  to  determine  the  average  pressure 
drop  through  the  cooling  jackets. 


Figure  19.  PDE  Cooling  Water  Supply  and  Discharge  Manifolds 


A  separate  system  provided  cooling  to  each  of  the  four  fuel  arms.  Water  was 
supplied  by  RPL  potable  water  and  the  piping  wrapped  around  the  fuel  arms.  The  piping 
was  then  enclosed  by  thermal  paste.  The  thermal  paste  can  be  seen  on  the  fuel  arms  in 
Figure  1 1 . 
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D.  INSTRUMENTATION 


1.  Thermocouples 

Type  K  thermocouples  from  Omega  Engineering,  Inc.  were  used  to  measure 
temperatures  at  various  locations.  Each  K  type  thennocouple  is  rated  up  to  1335  degrees 
C  with  a  response  time  of  0.55  seconds.  Each  thermocouple  was  wired  to  the  data 
acquisition  system  into  National  Instrument’s  TC-2095  rack-mounted  analog  breakout 
accessory.  Table  1  lists  all  the  thermocouple  temperatures  recorded  by  the  data 
acquisition  system. 


Table  1.  Thermocouple  Connections 


Thermocouple 

# 

Device 

National  Instrument 

TC-2095 

1 

Vitiator 

CH  00 

2 

Engine  Inlet 

CHOI 

3 

Cooling  Water  Inlet 

CH  02 

4 

Cooling  Water  Discharge  1 

CH  03 

5 

Cooling  Water  Discharge  2 

CH  04 

6 

Cooling  Water  Discharge  3 

CH  05 

7 

Cooling  Water  Discharge  4 

CH  06 

8 

Cooling  Water  Discharge  5 

CH  07 

9 

Cooling  Water  Discharge  6 

CH  08 

10 

Cooling  Water  Discharge  7 

CH  09 

11 

Cooling  Water  Discharge  8 

CH  10 

12 

Cooling  Water  Nozzle  Discharge 

CH  11 
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2. 


Pressure  Transducers 


Several  pressure  transducers  and  one  differential  pressure  cell  were  installed  to 
monitor  and  record  operating  of  the  PDE  and  cooling  water  systems.  The  pressure 
transducers  used  were  fabricated  by  Omega  Engineering,  Inc  Model  PX613.  Each 
transducer  is  ideal  for  gas  or  liquid  with  an  accuracy  of  0.4%.  Also,  one  Honeywell 
differential  pressure  cell,  Model  FDW,  was  connected  to  measure  the  differential 
pressure  across  the  cooling  water  manifolds.  All  pressure  cells  were  connected  to 
National  Instruments’  BNC  2090A  rack-mounted  analog  breakout  accessory.  Table  2 
lists  the  systems  measured  and  recorded  pressures  by  the  data  acquisition  system. 


Table  2.  Pressure  Transducer  Connections 


Pressure 

Transducer 

Range 

(GPa) 

Device 

National  Instrument 

BNC-2090A 

1 

0-13.8 

Main  Air  Supply 

AI7 

2 

0-13.8 

Ethylene  Supply 

AI8 

3 

0-13.8 

Hydrogen  Supply 

AI9 

4 

0-1.38 

Cooling  Water 

AI10 

5 

69000  Pa  differential 

Cooling  Water  D/P 

All  1 

6 

0-13.8 

Engine  Inlet 

AI12 

7 

0-1.38 

Shop  Air  Supply 

All  3 

3.  Kistler  Pressure  Sensors 

Two  Kistler  603B1  piezoelectric  pressure  transducers  were  inserted  into  the 
Kistler  228P  cooling  jackets.  Each  unit  was  installed  into  a  copper  spacer  placed 
between  to  combustor  sections.  The  copper  section  was  designed  by  Nichols  [8]  to  allow 
penetration  of  the  Kistler  probe  close  to  the  combustion  chamber  inner  wall  just  aft  of  the 
last  set  of  swept  ramp  obstacles  installed.  These  probes  are  capable  of  high  frequency 

pressure  measurements  under  the  extremely  high  operating  temperatures  of  detonation. 
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Figure  20.  Kistler  Pressure  Sensor  and  Kistler  Cooling  Jacket  (From  [8]) 

Mechanical  stress  due  to  a  pressure  wave  passing  by  each  of  the  Kistler  probes 
electrically  charges  the  piezoelectric  crystals,  which  then  provide  a  charge  proportional  to 
the  force  experienced  by  the  crystal.  The  charge  is  measured  in  picocoulombs  (pC)  and 
was  transmitted  to  a  Kistler  5010  multi-charge  amplifier  to  convert  and  amplify  the  signal 
to  a  proportional  voltage.  The  scale  of  the  amplifier  was  set  to  100  MU/volt  with  a 
sensitivity  of  0.380  pC/MU.  By  knowing  the  distance  between  the  two  probes,  the 
elapsed  time  measured  between  the  pressure  spikes  of  each  probe  can  determine  the 
speed  of  the  pressure  wave.  An  indication  that  detonation  occurred  comes  as  a  result  of  a 
wave  speed  in  excess  of  1,500  m/s. 


Figure  2 1 .  Kistler  Amplifiers 
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E.  DATA  ACQUISITION 


The  data  acquisition  was  recorded  by  two  National  Instrument  Graphical  User 
Interface  (GUI)  programs.  The  Kistler  high  speed  pressure  data  was  recorded  by  GUI 
pictured  in  Figure  22.  At  the  moment  the  PDE  was  ignited,  an  operator  would  depress 
the  “Start  Data  Recording”  and  three  seconds  of  high  speed  data  was  logged.  The  graph 
presented  in  the  Lab  View  program  displayed  the  results  that  allowed  immediate  analysis 
of  the  PDE  operation.  The  saved  data  was  later  analyzed  in  further  detail  to  ensure 
adequate  detonations  were  occurring. 


Figure  22.  National  Instruments  LabView  Data  Acquisition  System 


The  pressure  and  temperature  data  was  recorded  by  the  GUI  pictured  in  Figure  23. 
The  LabView  program  served  as  a  sequence  controller  for  the  operation  of  the  PDE  and 
also  as  data  acquisition.  Temperature  and  pressure  parameters  were  present  in  real  time 
on  the  screen  to  aid  the  operator  in  the  proper  operation  of  the  PDE,  and  also  were 
recorded  for  aiding  in  the  troubleshooting  of  the  PDC  operation  if  necessary.  The  logged 
pressure  and  temperature  data  was  analyzed  to  determine  corresponding  heat  transfer 
characteristics. 
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F. 


PDE  CONTROLLER  AND  OPERATING  PROCEDURE 


Two  National  Instruments  LabView  programs  served  as  the  controllers  for  the 
PDC.  Each  program  was  installed  on  separate  computers  and  controlled  by  two  operators 
in  the  RPL  control  room.  One  computer  was  connected  to  a  National  Instruments  PXIe- 
1062Q  controller  and  the  National  Instruments  PXI-1000B  controller  was  connected  to  a 
different  computer.  Together,  these  two  programs  set  the  operation  conditions  of  the 
PDE  and  remotely  controlled  the  gas  supply  valves  located  in  the  test  cell  although  the 
main  air  supply  and  cooling  water  flow  was  controlled  by  two  manual  switches  in  the 
control  room.  For  safety  considerations,  the  ignition  of  the  PDE  occur  manually  by  an 
operator  and  the  operator  was  capable  of  stopping  the  ignition  sequence.  For  further 
safety,  two  emergency  shutoff  buttons  are  located  in  the  control  room.  By  depressing  the 
shutoff  buttons,  fuel  injection  and  ignition  trigger  signals  were  disrupted  and  placed  the 
PDC  in  a  safe  condition  by  disabling  the  test  cell. 


Figure  23.  National  Instruments  LabView  Test  Cell  #2  PDE  Controller  1 
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The  GUI  pictured  in  Figure  24  sets  all  the  parameters  of  the  PDC  including 
frequency,  number  of  cycles,  fill  times  and  ignition  delay  times.  Once  engine  parameters 
were  set,  the  main  air  supply  was  initiated  followed  by  depressing  the  “Start  Sequence” 
button  on  the  GUI  in  Figure  23.  This  controller  operated  the  position  of  the  gas  valves, 
started  the  vitiator  sequence,  and  provided  the  countdown  to  manually  ignite  the  PDE.  A 
standard  operating  procedure  was  put  into  place  to  ensure  proper  valve  lineup  for  PDE 
startup,  proper  sequence  for  ignition,  and  to  provide  a  safe  condition  directly  after  a  run 
was  complete  as  well  as  for  a  complete  shutdown  of  the  PDE. 


Figure  24.  LabView  Interface  Controller  (From  [8]) 


G.  COOLING  WATER  MASS  FLOWRATE 

Testing  was  performed  on  the  cooling  water  system  to  calculate  the  mass  flow 

rates  of  the  cooling  jacket  sections  at  various  differential  pressures.  For  each  cooling 

jacket  section,  the  cooling  water  outlet  was  disconnected  from  the  cooling  water 
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discharge  manifold.  While  maintaining  a  constant  differential  pressure  across  the  cooling 
jacket,  cooling  water  was  collected  into  a  2000-milliliter  beaker  for  a  prescribed  amount 
of  time.  The  mass  of  the  water  in  the  beaker  was  measured  with  a  digital  scale.  This 
process  was  repeated  three  times,  and  the  mass  was  averaged  to  minimize  errors  due  to 
the  natural  inconsistencies  in  the  time  recording  evolution.  Each  cooling  jacket  was 
tested  in  this  manner  at  a  set  differential  pressure.  The  entire  process  would  then  be 
repeated  a  various  differential  pressures  by  adjusting  the  position  of  the  recirculation 
valve  in  the  cooling  water  system. 

Using  the  recorded  water  mass  for  each  calibration  point,  water  temperature,  and 
time,  mass  flow  rates  for  each  cooling  jacket  were  calculated  for  various  differential 
pressures.  It  was  observed  that  the  mass  flow  rates  for  each  individual  cooling  jacket 
were  varied  by  a  common  factor  between  the  two  differential  pressures. 

Since  the  mass  flow  rate  varied  by  the  same  amount  for  each  cooling  jacket,  only 
additional  measurements  were  taken  using  Cooling  Water  Jacket  #6  mass  flow  rate  at 
different  pressures.  After  plotting  the  mass  flow  rates  against  differential  pressure,  it  was 
noted  that  each  cooling  water  jacket  had  approximately  the  same  slope  for  a  given 
differential  pressure  change. 


Table  3.  Various  Cooling  Water  Mass  Flow  Rates 


Mass  Flow  Rate 

Differential  Pressure  (kPa) 

6.89 

13.8 

20.7 

27.6 

CW1 

(g/s) 

61.61 

72.32 

83.02 

93.73 

CW2 

(g/s) 

61.57 

72.28 

82.98 

93.69 

CW3 

(g/s) 

61.28 

71.99 

82.69 

93.40 

CW4 

(g/s) 

62.49 

73.20 

83.90 

94.61 

CW5 

(g/s) 

61.28 

71.99 

82.70 

93.40 

CW6 

(g/s) 

57.67 

68.37 

79.08 

89.79 

CW7 

(g/s) 

58.06 

68.77 

79.48 

90.18 

CW8 

(g/s) 

57.43 

68.13 

78.84 

89.55 

CWnozzle 

(g/s) 

48.57 

59.28 

69.98 

80.69 
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H.  PRESSURE  TRANSDUCER  CALIBRATION 


Calibration  of  the  pressure  transducers  was  necessary  to  ensure  proper  accuracy 
of  the  recorded  data  from  each  pressure  transducer.  Each  transducer  was  connected  to  a 
Heise  gauge  and  nitrogen  gas  bottle.  A  regulator  was  used  to  pressurize  each  transducer 
to  several  predetermined  set  points.  Once  pressure  was  stabilized  on  the  Heise  gauge,  the 
transducer  output  voltage  was  recorded.  The  recorded  data  was  used  to  develop  an 
equation  of  a  line  and  the  data  was  entered  into  the  Lab  View  program. 


Figure  25.  Calibration  of  Pressure  Transducers 


Table  4.  Pressure  Transducer  Calibrations 


Pressure 

Measurements 

Calibration  Data 

Cooling  Water 

P  =  50.18V  -  51.17 

Ethylene 

P  =  250.2V  -250.0 

Main  Air 

P  =  250.1V -250.0 

Engine  Inlet 

P  =  249.7V -251.6 
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IV.  EXPERIMENTAL  RESULTS 


The  PDC  was  operated  with  ethylene/air  mixtures  to  gather  heat  transfer 
characteristics  while  operating  at  different  frequencies  and  fuel  till  times.  By  keeping  the 
till  duration  constant  and  varying  frequency,  the  amount  of  time  that  purging  occurs 
relative  to  the  filling  of  fuel  and  air  was  altered  and  directly  affected  the  amount  of 
cooling  provided  by  the  purge  and  refresh  portion  of  the  cycle.  By  increasing  frequency, 
the  number  of  cycles  per  second  increased  the  aggregate  amount  of  heat  generated  by  the 
PDC.  Lastly,  by  changing  the  fuel  filling  times,  the  amount  of  fuel  available  for 
detonation  was  also  varied  thus  affecting  the  amount  of  heat  produced  for  each  fuel  cycle. 
By  testing  at  various  frequencies  and  fill  times,  the  variation  of  the  inner  combustor  wall 
was  observed  as  well  as  the  effective  heat  transfer  coefficient. 

Measuring  the  temperature  differential  across  each  cooling  jacket  provided  heat 
transfer  results  at  multiple  locations  along  the  combustor  axis.  Using  the  conduction  and 
convection  equations  from  Section  III,  the  cooling  water  temperatures  can  approximate 
the  combustor  wall  temperatures  along  its  entire  length. 

A.  TESTING  AT  20  HZ 

The  PDC  was  operated  at  a  frequency  of  20  Hz  with  a  stoichiometric  ethylene/air  mixture 
fill  time  of  20  milliseconds.  It  was  determined  the  cooling  water  outlet  temperatures 
reached  a  steady  state  condition  after  about  1 8  seconds  of  operation. 

In  Figure  26,  the  cooling  water  temperature  differentials  and  the  cooling  water 
mass  flow  rates  through  each  individual  cooling  jacket  was  used  to  determine  the  average 
heat  transfer  rate  in  Watts.  The  maximum  observed  heat  transfer  occurred  at  Cooling 
Water  Jackets  7  and  8.  These  two  cooling  water  jackets  are  positioned  just  aft  of  the  last 
swept  ramp  obstacle  in  the  combustor  tube.  The  position  of  the  last  ramp  has  been 
observed  during  previous  testing  to  be  the  location  where  detonation  initiation  occurs  and 
therefore  produces  high  temperature  and  pressures  from  the  detonation  event.  The 
maximum  heat  transfer  rate  is  observed  at  approximately  12  kilowatts  while  the 
minimum  is  around  9  kilowatts. 
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Heat  Transfer  Rate  vs  Time 


Figure  26.  Heat  Transfer  Rates  for  20  Hz  Tests 


Using  the  heat  transfer  rates  and  the  mean  cooling  water  temperature  for  each 
cooling  water  jacket,  the  liquid-side  wall  temperature  of  the  cooling  water  jacket  was 
calculated.  The  calculated  cooling  water  jacket  wall  temperature  was  then  used  with  the 
conduction  equations  to  calculate  the  inside  wall  temperature  of  the  combustor,  T  wall .  In 

Figure  27,  the  maximum  TgwalI  occurs  at  Cooling  Water  Jackets  7  and  8,  which 
corresponds  to  the  maximum  heat  transfer  rates.  The  maximum  combustor  wall 
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temperature  was  determined  to  be  slightly  less  than  130°  C.  The  cooling  water 
differential  pressure  was  set  to  approximately  20,700  Pa  during  the  20  Hz  tests.  This 
flow  rate  is  achieved  by  setting  the  recirculation  line  to  full  open  position  allowing  the 
minimum  flow  rate  without  adjusting  the  bah  valve  leading  to  Test  Cell  #2.  A  schematic 
of  the  Cooling  Water  System  is  located  in  Appendix  B.  The  amount  of  cooling  water 
mass  flow  provides  more  than  adequate  cooling  of  the  combustor  inner  wall. 


Combustor  Wall  Temperature  vs.  Time 


Figure  27.  Combustor  Wall  Temperatures  for  20  Hz  Tests 


Figure  28  displays  the  variations  in  combustor  wall  temperatures  along  the  length 
of  the  combustor  for  time  steps  of  0,  5,  10,  15  and  20  seconds.  The  shape  of  the 
temperature  profile  remains  nearly  constant  from  the  initial  start  to  steady  state 
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conditions.  An  increasing  temperature  is  noticed  from  near  the  front  of  the  combustor 
towards  the  end  as  DDT  occurs  with  a  maximum  just  prior  to  reaching  the  nozzle  inlet. 


B.  TESTING  AT  30  HZ 

Testing  at  30  Hz  with  the  same  fill  fraction  produces  50%  more  detonations  per 
second  as  compared  to  the  20  Hz  testing.  With  many  more  detonations  occurring,  a 
noticeable  increase  in  heat  transfer  rates  and  combustor  wall  temperature,  T  wall,  is 

observed.  It  is  assumed  that  the  increase  in  the  number  detonations  per  second  caused 
some  internal  components  in  the  combustor  to  reach  temperatures  above  the  auto  ignition 
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point  for  ethylene  gas/air  mixture.  Once  the  auto  ignition  temperatures  were  present, 
flame  holding  was  observed  in  the  PDC.  The  components  assumed  to  be  the  source  for 
auto  ignition  were  either  the  ramp  obstacles  down  the  combustion  chamber,  the  shroud 
that  surrounds  the  Transient  Plasma  Igniter  or  a  thermocouple  directly  in  the  flow  field 
near  the  TPI  unit. 

The  four  sets  of  swept  ramp  obstacles  were  located  at  Cooling  Water  Jacket 
sections  3-6.  The  ramps  are  held  in  place  by  a  retaining  bolt  place  through  the  flange  of 
the  cooling  jacket.  It  was  assumed  that  a  small  uniform  layer  of  silicone  existed  between 
the  ramp  and  the  combustor  wall  due  to  the  assembly  technique,  therefore  minimizing  the 
amount  of  heat  transfer  possible  from  the  ramps  to  the  cooling  jacket.  The  thermocouple 
was  also  removed  from  the  flow  field  to  eliminate  its  potential  as  the  source  of  the  auto 
ignition.  The  shroud  located  at  the  front  of  the  combustor  is  located  along  the  centerline 
of  the  ignition  region.  The  length  of  the  shroud  was  cut  to  shorter  lengths  to  attempt  to 
provide  improved  cooling  and  minimize  the  chance  of  auto  ignition.  As  the  length  of  the 
shroud  was  shortened,  the  length  of  time  prior  to  auto  ignition  occurring  increased  from 
about  5  seconds  to  10  seconds.  The  source  of  auto  ignition  could  not  be  isolated  to  either 
the  TPI  shroud  or  the  swept  ramp  obstacles  so  tests  at  30  Hz  was  limited  to 
approximately  10  seconds. 

Initial  30  Hz  testing  occurred  with  a  fuel/air  fill  time  of  18  milliseconds,  which 
corresponds  to  fully  filling  the  combustion  chamber.  In  Figure  29,  the  maximum  heat 
transfer  rates  occurred  at  Cooling  Water  Jackets  7  and  8,  which  is  consistent  with  the  20 
Hz  testing.  The  maximum  observed  heat  transfer  rates  reached  approximately  15 
kilowatts.  Unfortunately,  the  heat  transfer  rates  were  not  able  to  reach  steady  state 
conditions  prior  to  the  occurrence  of  auto  ignition.  The  rates  after  10  seconds  of 
operation  are  greater  than  the  steady  state  conditions  reached  at  20  Hz  testing.  The  rate 
of  change  of  the  heat  transfer  rates  at  10  seconds  appears  to  be  decreasing,  which  means 
it  is  attempting  to  start  to  level  off  to  steady  state  conditions.  It  is  assumed  that  steady 
state  conditions  would  be  reached  in  less  than  20  seconds. 


39 


Heat  Transfer  Rate  vs  Time 


Figure  29.  Heat  Transfer  Rates  for  30  Hz  Testing 


In  Figure  30,  the  heat  transfer  rate  profiles  along  the  length  of  combustion 
chamber  for  30  Hz  tests  have  the  same  shape  as  seen  in  the  20  Hz  tests.  The  profiles 
shown  are  in  two  second  intervals  from  0-10  seconds  with  steadily  increasing  heat 
transfer  rates  along  each  time  step.  The  front  end  of  the  combustor  is  producing  about  1 1 
kilowatts  while  the  aft  end  where  the  DDT  event  occurs,  experiences  rates  as  high  as  15 
kilowatts  prior  to  reaching  a  flame  holding  condition. 
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Figure  30.  Heat  Transfer  Rate  Profiles  for  30  Hz  Tests 

The  fuel  fill  times  were  varied  from  14-18  milliseconds  for  30  Hz  operation  to 
observe  the  changes  in  heat  transfer.  At  18  ms,  fuel  completely  fills  the  combustion 
chamber  and  results  in  the  greatest  heat  transfer  rates  and  maximum  combustor  inner  wall 
temperatures.  Figure  31  displays  a  comparison  of  fill  times  of  14,  16,  and  18  ms  at 
Cooling  Water  Jacket  7.  All  three  fill  times  produce  combustor  wall  temperatures  higher 
than  the  20  Hz  tests  and  it  is  obvious  that  the  temperatures  are  still  increasing  and  have 
not  reach  steady  state  in  all  three  conditions.  Even  though  steady  state  conditions  were 
not  achieved,  it  appears  as  if  the  inner  combustor  wall  temperature  would  remain  less 
than  180  degrees  C.  The  outcome  confirms  that  a  partial  fill  of  fuel  in  the  combustion 

chamber  results  in  a  lower  amount  of  heat  transfer  at  the  aft  end  of  the  combustor. 
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Combustor  Wall  Temperature  vs.  Time 


Figure  31.  Combustor  Wall  Temperatures  at  CW7  for  Fuel  Fill  Times  of  14-18  ms 


C.  COMPARISONS  OF  20HZ  AND  30  HZ  TESTING 

The  results  from  20  Hz  and  30  Hz  testing  were  compared  since  steady  state 
conditions  are  reached  during  20  Hz  testing  and  not  for  30  Hz.  Comparisons  were  made 
with  a  fuel  fill  time  of  18  ms,  which  corresponds  to  a  full  fill  of  the  combustion  chamber. 

In  Figure  33,  the  temperatures  of  20  Hz  and  30  Hz  testing  at  Cooling  Water  Jacket 
7  are  compared.  The  20  Hz  testing  occurred  for  20  seconds  while  the  30  Hz  testing  was 
stopped  after  10  seconds  of  operation.  The  30  Hz  temperature  at  10  seconds  reaches 
approximately  150  °C,  which  is  about  30%  greater  than  the  20  Hz  temperature  after  10 
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seconds  of  operation.  It  appears  as  if  the  30  Hz  testing  were  to  reach  steady  state 
conditions  then  the  inner  combustor  wall  temperature  would  not  exceed  200°  C.  At 
these  temperatures,  the  combustor  wall  would  be  cool  enough  for  extended  operation. 


Time  (s) 


Figure  32.  Combustor  Wall  Temperature  Comparisons  at  20  Hz  and  30  Hz 


The  temperature  profiles  along  the  length  of  the  combustion  chamber  were 
compared  at  three  second  intervals  for  20  Hz  and  30  Hz  testing.  In  Figure  34,  the  dotted 
lined  lines  represent  the  30  Hz  testing  temperatures  while  the  solid  lines  are  for  20  Hz 
temperatures.  The  shape  of  the  profiles  remains  similar  for  each  time  step.  The  graph 
displays  an  increasing  separation  of  temperatures  from  20  Hz  to  30  Hz  as  the  run 
increases.  The  30  Hz  temperatures  are  at  least  40%  greater  than  the  20  Hz  temperatures. 
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A  test  run  was  conducted  without  any  ignition  and  only  used  the  vitiator  to  heat 
the  incoming  air  into  the  combustor  to  allow  the  initial  equilibrium  conditions  to  exits. 
On  Figure  34,  the  combustor  wall  temperature  from  only  the  vitiator  nearly  duplicates  the 
initial  conditions  for  the  start  of  the  20  Hz  and  30  Hz  testing. 


Figure  33.  Temperature  Profile  Comparisons  at  20  Hz  and  30  Hz 


In  Figure  35,  the  heat  transfer  rates  at  Cooling  Water  Jacket  7  are  compared  for  20 
Hz  and  30  Hz  testing.  The  20  Hz  heat  transfer  reached  a  maximum  at  about  12  kilowatts 
at  an  equilibrium  condition.  The  30  Hz  data  was  plotted  for  10  seconds  of  operation  and 
then  a  projected  profile  was  continued  to  beyond  20  seconds  along  the  same  shape.  This 
projected  profile  produces  a  maximum  heat  transfer  rate  of  18  kilowatts,  which  is  50% 


44 


greater  than  the  maximum  for  20  Hz  testing.  This  increase  in  heat  transfer  corresponds  to 
a  50%  increase  in  the  number  of  detonations  that  occurs  per  second. 


Heat  Transfer  Rate  vs  Time 


Figure  34.  Combustor  Wall  Temperature  Comparisons  at  CW7  for  20  Hz  and  30  Hz 


D.  COMPARIONS  OF  30  HZ  AND  40  HZ  TESTING 

Several  tests  were  conducted  at  40  Hz  and  the  results  were  compared  to  30  Hz 
testing.  At  a  frequency  of  40  Hz,  33%  more  detonations  occur  than  during  a  frequency  of 
30  Hz  resulting  in  an  increased  heat  transfer  rate.  The  increased  frequency  resulted  in 
heat  transfer  rates  as  high  as  16  kilowatts  prior  to  reaching  10  seconds  of  testing  as  well 
as  higher  inner  combustor  wall  temperatures.  Figure  32  displays  a  comparison  of 
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combustor  inner  wall  temperatures  at  CW  7  for  30  Hz  and  40  Hz  testing.  It  appears  that 
the  40  Hz  temperatures  are  trending  about  15%  greater  than  30  Hz  temperatures. 


Time  (s) 


Figure  35.  Combustor  Wall  Temperatures  at  CW7  for  30  and  40  Hz  Tests 


E.  COMPUTATIONAL  FLUID  DYNAMICS  MODELING 

Computational  fluid  dynamic  modeling  was  performed  using  ANSYS  CFX 
software  to  reproduce  the  temperature  profiles  that  occur  during  detonation  with  an 
isothermal  combustor  wall.  The  CFD  was  setup  for  a  0.05  meter  diameter  combustion 
chamber  without  swept  ramps  using  stoichiometric  conditions  for  ethylene  gas  and  air 
with  an  initial  temperature  of  530  K  and  pressure  of  506,625  Pa  (5  atm).  Figure  36 
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shows  four  different  time  steps  during  the  detonation  event  in  the  combustor.  The  first 
image  is  at  250  ms  where  the  detonation  wave  has  just  passed  point  B,  about  three- 
quarters  the  length  of  the  combustor.  At  point  B,  the  temperatures  reached  were  as  high 
as  3,000  K  and  close  to  the  Chapman- Jouguet  conditions.  However,  the  temperature  does 
not  stay  this  high.  As  the  detonation  wave  passes,  the  temperature  quickly  drops  over 
600  K  in  about  0.5  ms  as  seen  in  the  second  image.  In  the  second  image,  the  detonation 
has  reached  the  nozzle  and  some  shock  wave  are  reflected  and  heading  towards  the  front 
of  the  combustor  during  the  phase  known  as  ring  down.  The  shock  wave  essentially 
“resets”  the  boundary  layer  conditions.  As  more  time  passes,  the  boundary  begins  to 
grow  but  is  still  disturbed  at  1.5  ms  after  the  initial  detonation  wave  passes.  In  the  last 
image,  the  shock  waves  are  further  dissipating  and  the  combustor  is  nearly  back  to 
refresh  conditions.  These  conditions  are  highly  transient  and  make  determining  a  heat 
transfer  coefficient  very  difficult. 


47 


L 


Ethylene  Air  Detonation  (Phi  =  1.0) 

28Feb2011 

Pinitial  =  506625  Pa 

Tinitial  =  530  K 

Twall  =  480  K  (Isothermal) 


0.250  ms 


B  A 


T:  600  800  1000  1200  1400  1600  1800  2000  2200  2400  2600  2800  3000 


0.750  ms 


■ - ! 


1.250  ms 


J - L 


Figure  36.  CFD  Temperature  Profiles  for  Detonation 


Figure  37  displays  a  comparison  of  the  temperature  and  pressure  variations  along 
the  combustor  wall  at  points  A,  B,  and  C.  At  point  A,  the  refresh  conditions  deliver  low 
temperatures  and  pressures  in  the  flow  field.  Immediately  behind  the  detonation  wave  at 
point  B  (red  lines),  a  rapid  jump  in  temperature  and  pressure  is  observed  right  up  to  the 
combustor  wall  with  a  very  small  boundary  layer  thickness.  Further  back  at  point  C, 
pressure  and  temperature  in  the  bulk  flow  have  substantially  reduced. 
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Ethylene  Air  Detonation  (Phi  =  1 .0) 
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Figure  37.  Temperature  and  Pressure  Profile  in  a  Detonation  Wave  (0.250  ms) 


Figures  38  and  39  show  the  boundary  layer  temperature  profiles  at  various  stages 
during  the  ring  down  event  after  detonation  has  occurred.  The  average  bulk  flow 
temperature  was  determined  to  be  approximately  2,400-2,500  degrees  K.  This  provides 
about  a  2,000  degree  K  differential  temperature  from  the  bulk  flow  conditions  to  the 
combustor  inner  wall.  This  increase  in  thermal  and  pressure  conditions  from  the 
detonation  are  the  main  driving  force  of  the  heat  transfer  across  the  combustor  wall  that 
occurs  for  very  short  time  frame  of  3—4  milliseconds.  After  the  ring  down  event  occurs, 
the  refresh  portion  of  the  cycle  begins  and  the  combustion  chamber  is  filled  with  more  air 
and  eventually  more  fuel. 
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Ethylene  Air  Detonation  (Phi  =1.0) 
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Figure  38.  Temperature  Profile  in  a  Detonation  Wave  (1.250  ms) 
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Ethylene  Air  Detonation  (Phi  =  1.0) 
28Feb2011 

Pinitial  =  506625  Pa 

Tinitial  =  530  K 

Twall  =  480  K  (Isothermal) 

Figure  39.  Temperature  Profile  in  a  Detonation  Wave  (1.750  ms) 

In  order  to  determine  an  effective  heat  transfer  coefficient,  the  detonation  event 
conditions  have  to  be  accounted  for  as  well  as  the  refresh  conditions.  During  30  FIz 
operation,  one  cycle  occurs  every  33.3  ms  with  the  detonation  event  lasting  for 
approximately  10%  of  the  cycle  time.  Refresh  conditions  can  be  assumed  for  the 
remaining  cycle  time.  An  average  heat  transfer  coefficient  was  determined  by  using  the 
conditions  from  the  CFD  and  the  experimentally  measured  data  from  the  refresh 
conditions  using  the  vitiator.  For  20  FIz  operation,  one  cycle  lasts  for  50  ms  with  3^1  ms 
for  the  detonation  event  and  the  remaining  time  involves  purging/filling  the  combustor 
with  fuel  or  air.  With  the  maximum  heat  transfer  rate  equivalent  to  12,000  Watts  for  20 

FIz  steady  state  testing,  one  cycle  transfers  600  J  of  energy  to  a  single  cooling  jacket. 
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During  the  refresh  conditions,  138  J  is  delivered  based  on  experimentally  determined 
data.  The  detonation  event  must  deliver  462  J  in  4  ms,  resulting  in  a  heat  transfer  rate  of 
115,500  W.  Using  the  surface  area  that  the  heat  transfer  occurs  across  and  the 
differential  temperature  determined  from  the  CFD,  the  heat  transfer  coefficient  during 
detonation  would  be  approximately  3,166  W/nrK  .  This  would  result  in  an  average  in 
the  883  W/m2K  .  While  assuming  the  maximum  heat  transfer  rate  of  18,000  W  for  30  Hz 
testing,  the  same  calculation  can  be  repeated  adjusting  the  cycle  time  and  results  in  an 
average  heat  transfer  coefficient  of  1,035  W/m2K  . 

It  was  calculated  that  each  PDC  cycle  produces  almost  300  kilowatts  using 
ethylene  at  a  pressure  of  250,000  Pa  and  a  temperature  of  400  degrees  K,  which  is  similar 
to  operating  conditions  conducted  during  testing.  Using  this  information  determined  that 
about  1.6%  of  the  energy  released  during  combustion  produced  was  transferred  to  the 
cooling  water  system.  This  transferred  energy  could  possibly  be  reused  in  a  regenerative 
system  to  preheat  the  incoming  air  or  fuel. 

F.  UNCERTAINTY  ANALYSIS 

An  analysis  was  made  to  determine  the  total  error  associated  with  equipment 
accuracies  and  human  error.  The  type  K  thermocouples  used  to  measure  cooling  water 
temperatures  had  an  accuracy  range  of  ±  2.2  °  C  for  temperature  below  293  °  C.  A  type  K 
thermometer  measured  the  inlet  and  outlet  temperature,  therefore  doubling  the  2.2  °C 
possible  inaccuracy.  A  0.1-second  response  delay  time  was  used  for  calculating  the  error 
in  the  mass  flow  rate  measurements.  A  ±  690  Pa  inaccuracy  was  assumed  for  human 
error  of  observing  the  cooling  water  system  differential  pressure.  These  combined 
inaccuracies  and  errors  resulted  in  an  approximate  error  of  16.4%,  which  results  in  an 
approximate  heat  transfer  rate  of  11,500  W  ±  1,900  W.  The  type  K  thermocouples 
accounted  for  the  majority  of  total  error  calculated.  For  future  studies,  a  type  T 
thermocouple  should  be  used  with  an  accuracy  of  ±  1°C.  For  an  even  better  result,  a 
Resistance  Temperature  Detector  (RTD)  would  provide  an  accuracy  of  less  than  ±  0.2  °  C 
while  operating  below  100  °C. 
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The  accuracy  of  the  calculated  combustor  inner  wall  temperature  was  also 
determined.  The  Nusselt  number  was  determined  from  Equation  11,  the  Dittus-Boelter 
equation,  which  may  result  in  errors  as  large  as  25%  [7],  A  generally  more  complex 
equation  attributed  to  Petukhov  may  reduce  the  error  to  less  than  10%  is  shown  in 
Equation  18,  which  uses  the  friction  factor  f,  from  the  Moody  diagram  [7], 


NuD 


// 


Re^  Pr 


1.07  +  12.7 


// 


(Pr73-1) 


(18)  (From  [7]) 


The  variation  in  thermal  conductivity  for  4340  steel  was  considered  negligible  for 
the  given  operating  conditions.  An  additional  ±2.2°  C  was  added  to  the  mean 
temperature  in  Equation  13  for  the  K  type  thermocouple.  The  combined  error  associated 
with  the  heat  transfer  rate,  Nusselt  number,  and  thermocouple  resulted  in  a  ±  12.0%  total 
error  in  the  calculated  inner  combustor  wall  temperature  Tg  wall  at  steady  state  20  Hz 


operations.  The  calculated  T  a  of  182°C  may  vary  as  much  as  ±21.9°C.  This  error 
could  be  reduced  by  using  Equation  18  and  a  more  accurate  thennocouple. 


The  effects  of  natural  convection  and  radiation  were  calculated  for  the  combustor 
operating  in  a  23  °C  environment  assuming  the  exterior  cooling  jacket  wall  reach 
equilibrium  temperature  with  the  temperature  of  the  cooling  water.  With  assuming  an 
emissivity  of  0.85,  each  cooling  jacket  section  would  have  an  addition  10-20  W,  which  is 
negligible  considering  the  orders  of  magnitude  greater  heat  transfer  occurs  between  the 
combustor  wall  and  cooling  jacket. 


53 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


54 


V.  CONCLUSIONS  AND  FUTURE  WORK 


A.  CONCLUSIONS 

This  thesis  explored  the  heat  transfer  characteristic  of  a  water-cooled  PDC 
through  experimental  testing  and  computational  analysis.  Operational  frequencies  from 
20-40  Hz  were  considered  and  found  to  be  directly  proportional  to  the  heat  transfer  rates 
to  the  combustor  wall.  By  varying  the  fill  time,  or  fill  fraction,  of  the  fuel  from  14-18 
milliseconds,  the  amount  of  energy  transferred  along  the  axis  was  found  to  noticeably 
vary  along  the  combustor. 

The  cooling  water  system  provided  adequate  cooling  to  the  combustor  wall  to 
allow  steady  state  conditions  to  be  achieved.  The  amount  of  cooling  necessary  to  cool  to 
the  combustor  wall  was  found  to  be  low  enough  such  that  cooling  could  be  achieved  by 
using  the  air  and/or  fuel  as  a  cooling  medium. 

The  results  from  experimental  tests  revealed  that  the  maximum  heat  transfer 
occurs  aft  of  the  swept  ramp  obstacles.  The  region  of  increased  heat  transfer  comes  as  a 
result  from  the  DDT  event  occurring  just  after  the  last  obstacle.  This  region  requires  the 
most  amount  of  cooling  and  should  be  considered  in  future  cooling  water  system  designs. 

The  transient  conditions  during  each  detonation  cycle  created  difficulties  in 
obtaining  an  instantaneous  heat  transfer  rate  as  well  as  an  instantaneous  heat  transfer 
coefficient.  By  examining  the  refresh/purge  condition  separate  from  the  detonation 
event,  an  average  effective  heat  transfer  coefficient  for  each  portion  of  the  cycle  can  be 
approximated  providing  valuable  information  towards  continuing  PDE  research. 
Computer  simulations  were  performed  to  provide  understanding  of  heat  transfer 
conditions  during  the  detonation  wave  cycle  and  allowed  for  the  investigation  of  the  flow 
field  characteristics  along  the  combustor  wall  as  compared  to  the  flow  field  conditions 
during  detonation  and  ring  down  conditions.  By  combining  the  experimental  research 
data  and  the  computer  simulations,  an  average  heat  transfer  coefficient  could  be 
approximated  for  different  operation  frequencies. 


55 


B.  FUTURE  WORK 

A  redesign  of  the  cooled  combustor  chamber  is  currently  in  progress  and  will 
result  in  a  cast  mold.  A  water-cooled  combustor  with  casted  swept  ramps  in  the 
combustion  chamber  that  are  cooled  as  well  maximizes  the  amount  cooling  to  the  ramps 
to  help  reduce  ramp  temperatures.  By  reducing  ramp  temperatures,  the  auto  ignition  of 
the  reactants  can  possibly  be  eliminated  as  being  caused  by  the  obstacles  in  the  flow  field. 

Further  research  can  be  conducted  using  different  fuels  such  as  liquid  JP-10, 
which  has  a  higher  flashpoint  than  ethylene/air  mixtures.  The  higher  flashpoint  could 
surpass  the  possible  flame-holding  event,  and  allow  for  longer  PDE  operation  at  higher 
frequencies. 

A  redesign  of  the  cooling  water  system  as  a  regenerative  heat  exchanger  using  air 
or  fuel  would  be  desirable.  A  regenerative  heat  exchange  would  allow  for  cooling  of  the 
combustor  wall  while  preheating  the  incoming  air  or  fuel  to  a  higher  initial  operating 
temperature.  For  an  air  system,  a  design  with  a  large  number  of  heat  exchanger  fins 
might  be  necessary  to  provide  adequate  heat  transfer  could  be  detennined  by  heat 
exchanger  effectiveness  calculations.  A  heat  exchanger  material  selection  study  would 
provide  relevant  data  to  which  material  can  withstand  the  high  frequency  mechanical 
stresses  and  the  associated  thermal  stresses  from  detonation  while  providing  adequate 
heat  transfer  properties. 
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APPENDIX  A:  PDE  STANDARD  OPERATING  PROCEDURES 


Standard  Operating  Procedures 

Test  Cell  #2 

Modification  Date  (27  January  2011) 

RUN  Setup  Procedures 

1 .  Lab  Personnel  -  NOTIFY  OF  IMPENDING  TEST 

2.  Gate -LOCK 

3.  Red  Emergency  Stop  Buttons  (x2)  -  IN 

4.  5V  Power  Supply  -  OFF 

5.  Verify  24V  Power  Supply  (Control  Room)  -  OFF 

6.  Warning  Lights  -  ON 

7.  Air  Bank  Pressure  -  CHECK  >1500  psi 

8.  Run  Sheet  -  COMPLETE 

a)  Required  pressures  -  NOTE 

9.  On  TC#3  Computer  (32-bit) 

a)  “TC2  PDE  Control  24  Jan”  -  OPEN  &  RUN 

i.  Change  data  file  name,  right  click  data  tile,  select  “Data 
Operations”,  select  “Make  Current  File  Default”,  File  -  SAVE 

b)  “PDE  High  Speed  30  July”  -  OPEN  &  RUN 

i.  Change  data  tile  name,  right  click  data  tile,  select  “Data 

Operations”,  select  “Make  Current  File  Default”,  File  -  SAVE 

10.  On  TC#2  Computer  (32-bit) 

a)  “National  Instruments  Lab  View”  -  OPEN 

b)  “Test  Cell  #2.1vprojecf  ’  -  OPEN 

c)  Maximize  tree  by  clicking  +  symbol 

d)  “Test  Cell  #2  with  Brady  Revamp”  -  OPEN  &  RUN 

e)  Enter  Values  from  Run  Sheet 

f)  “Set  Engine  Parameters”  -  SELECT 

1 1 .  BNC  Cabinet  Power  Strip  -  ON 

12.  BNC  Box  (on  top  of  cabinet)  -  ON 

a)  CH.  A  (0.00007  /  0.0)  -  VERIFY  (set  with  TC#2  computer) 

b)  CH.  B  (0.00005  /  0.0002 1)  -  VERIFY  (set  with  TC#2  computer) 

Outside 

13.  Jamesbury  Valve  -  OPEN 

14.  24  Volt  Power  Supply  (in  TC#0)  -  ON 

15.  Node  4  Ball  Valve  (in  TC#1)  -  OPEN 

16.  H2  Six  Pack  3  Bottle  Isolation  Valves  -  OPEN 

17.  H2  Six  Pack  3  Isolation  Valve-OPEN  &CHECK  PRESSURE 

18.  DAQ  Power  (in  TC#3)  -  ON 

19.  Cooling  Water  System 
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a)  Test  Cell  #3  Knife  Switch  -  ON 

b)  Knife  Switch  Breaker  Handle  -  ON 

c)  Water  Tank  -  CHECK  (full  and  clean) 

d)  Water  Tank  Isolation  Valve  -  OPEN 

e)  Test  Cell  #2  Ball  Valve  -  OPEN  (ensure  TC#3  valve  closed) 

20.  At  Overhead  Boxes  (in  TC#2) 

a)  Power  Supply  -  ON  (170  volts) 

b)  TPI  -  ON 

2 1 .  Vitiator  Spark  Plug  -  DISCONNECT 

22.  Main  Air  (yellow  handle)  -  CLOSE 

23.  Liquid  Injector  Cooling  Water  Valve  -  OPEN 

24.  Shop  Air  (red  handle)  -  OPEN  (can  verify  with  blue  handle) 

25.  Node  4  Isolation  Valve  -  OPEN 

26.  Kistler  Amps 

a)  Power  Switch  -  ON 

b)  Ensure  Proper  Setting  -  0.380mV 

c)  Function  Button  -  OPERATE 

27.  Transducer  TESCOM  Power  -  ON 

28.  Gas  pressure  on  Node  22  (N2)  to  appropriate  level  to  prevent  excessive  venting  - 
SET 

29.  Pressurized  Gases 

a)  Ethylene  Ball  Valve  -  OPEN 

i.  Check  C2H4  pressure  in  accumulator  and  note  if  sufficient.  If  NOT 
sufficient  perform  accumulator  fill  procedures 

b)  Ethylene  Ball  Valve  -  CLOSE 

c)  H2  -  OPEN 

d)  H2  Torch -OPEN 

e)  N2  Tank -OPEN 

30.  Shop  Air  Tank  (closet)  -  CHECK  (95-120  psi) 

Inside 

3 1 .  Set  Gas  Pressures  (in  control  room) 

a)  Node  1 ;  Main  Air 

b)  Node  4;  High  Pressure  Air 

c)  Node  20;  Vitiator  H2 

d)  Node  22;  C2H4  controlled  with  N2 

i.  Set  to  2x  C2H4  pressure 

32.  24  volt  DC -ON 

33.  BNC  Box -RUN 

Outside 

34.  Main  Air  (yellow  handle)  -  OPEN 

35.  Vitiator  Spark  Plug  -  CONNECT 
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*TEST  CELL  DANGER  CONDITION 


Run  Profile 


1.  Personnel  -  HEAD  COUNT 

2.  Lab  View  Programs  -  MODIFY  FILE  NAME  AS  NECESSARY  &  RUN 

3.  Golf  Course  -  CLEAR 

4.  Siren  -  ON 

5.  Emergency  Stop  Buttons  (x2)  -  OUT 

6.  5V  Power  Supply  -  ON 

7.  Main  Air  -  ON 

8.  Test  Cell  2  PDE  Controller  -  START  SEQUENCE 

9.  Cooling  Water  Pump  Switch  -  ON  (After  3-way  valve  operates) 

10.  Countdown  (Beginning  with  4  Second  Count  From  After  Time  Set  for  Vitiator) 

1 1 .  Bottom  BNC  Controller  -  START 

12.  PDE  High  Speed  Data  Recording  -  START 

After  Run 


1 .  Main  Air  -  OFF 

2.  Cooling  Water  Pump  Switch  -  OFF 

3.  Emergency  Stop  Buttons  (x2)  -  IN 

4.  5V  Power  Supply  -  OFF 

5.  Siren -OFF 

Run  Shutdown  Procedure 


1 .  Emergency  Stop  Buttons  (x2)  -  VERIFY  IN 

2.  5V  Power  Supply  -  VERIFY  OFF 

3.  Set  Gas  Pressures 

a.  Node  1  -  ZERO 

b.  Node  4 -ZERO 

c.  Node  20 -ZERO 

d.  Node  22  -  MAINTAIN  CURRENT  VALUE  (consider  minor  reduction) 

4.  BNC  Box -OFF 

5.  BNC  Cabinet  Power  Strip  -  OFF 

6.  24  volt  DC  -  OFF  (check  with  other  test  cells  prior) 

7.  Jamesbury  Valve  -  CLOSE 

8.  24  Volt  Power  Supply  (in  TC#0)  -  ON 

9.  Node  4  Ball  Valve  (in  TC#1)  -  CLOSE 

10.  Vitiator  Spark  Plug  -  DISCONNECT 

1 1 .  Main  Air  (yellow  handle)  -  CLOSE 

12.  Kistler  Amplifiers  -  OFF 

13.  At  Overhead  Boxes  (in  TC#2) 

a.  TPI-OFF 
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b.  Power  Supply  -  OFF 

14.  Water  Valve  -  CLOSE 

15.  Shop  Air  (red  handle)  -  CLOSE 

16.  Bleed  Shop  Air  (blue  handle)  -  OPEN  then  CLOSE 

17.  Node  4  Isolation  Valve  -  CLOSE 

18.  Transducer  TESCOM  Power  -  OFF 

19.  Pressurized  Tanks 

a.  H2- CLOSE 

b.  H2  Torch  -  CLOSE 

c.  N2- CLOSE 

20.  Cooling  Water  Pump 

a.  Test  Cell  #2  Ball  Valve  -  CLOSED 

b.  Water  Tank  Isolation  Valve  -  CLOSED 

c.  Knife  Switch  Breaker  Handle  -  OFF 

d.  Test  Cell  #3  Knife  Switch  -  OFF 

2 1 .  DAQ  Power  (in  TC#3)  -  OFF 

22.  H2  Six  Pack  -  CLOSE  &  RECORD  PRESSURES 

23.  Warning  Lights  -  OFF 
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APPENDIX  B:  COOLING  WATER  SYSTEM 


D/P 

Cell 


Figure  40.  Cooling  Water  System  Schematic 
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APPENDIX  C:  COMBUSTOR  ASSEMBLY  DRAWINGS 
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Figure  41.  Combustor  Section  -  Isometric  View  (From  [8]) 
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Figure  42.  Combustor  Sections  -  Plan  View  (From  [8]) 
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Figure  43.  Combustor  Sections  -  Inner  Tube  (From  [8]) 
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Figure  44.  Cooling  Nozzle  -  Inner  Tube  (From  [8]) 
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Figure  45.  Cooling  Nozzle  -  Assembly  (From  [8]) 
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